Persistence in the zero-temperature dynamics of the Q-states 
Potts model on undirected-directed Barabasi- Albert networks 
and Erdos-Renyi random graphs * 

F. P. Fernandes and F.W.S. Lima 

Departamento de Fisica, Universidade Federal do Piaui, 64049-550, Teresina 
- PI, Brazil 

e-mail: fwslima@gmail.com 

•k This paper is dedicated to Dietrich Stauffer for his great contribution 
to physics. 

Abstract: The zero-temperature Glauber dynamics is used to investigate the 
persistence probabihty P{t) in the Potts model with Q = 3,4, 5, 7, 9, 12, 24, 64, 128, 
256, 512, 1024, 4096, 16384,..., 2^0 states on directed and undirected Barabasi- Albert 
networks and Erdos-Renyi random graphs. In this model it is found that P{t) 
decays exponentially to zero in short times for directed and undirected Erdos- 
Renyi random graphs. For directed and undirected Barabasi-Albert networks, in 
contrast it decays exponentially to a constant value for long times, i.e, P(oo) is 
different from zero for all Q values (here studied) from Q = 3, 4, 5, 2^*^; this 
shows "blocking" for all these Q values. Except that for Q = 2^° in the undirected 
case P{t) tends exponentially to zero; this could be just a finite-size effect since in 
the other "blocking" cases you may have only a few unchanged spins. 
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Introduction 

Derrida et al. yij discovered that the fraction P{t) of spins which have 
not flipped until time t exhibits non-trivial dynamical behavior for the Potts 
model and zero-temperature dynamics. In one dimension, they showed that 
P{t) decays algebraically, with an exponent 6 varying continuously between 
and 1 as Q increases from 1 to oo [HE]- Derrida et al. [3] studied the 
problem in two dimensions by Monte Carlo simulations. However, their re- 
sult for high Q was inaccurate since the data showed some curvature and 
only effective exponents for 6 could be determined. Michael Hennecke |1] 
investigated the two-dimensional Potts model for Q = oo by a similar zero- 
temperature Monte Carlo method. Due to a much higher numerical effort. 
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its has been confirmed that asymptotically, 6 approaches unity if Q = oo. 
Michael Hennecke [S] presented also simulation results for finite Q in the 
order to obtain precise, and hopefully asymptotic, values for the exponents 
9 for general Q. He applied the Monte Carlo Method described in Ref. |3] to 
a triangular lattice with nearest-neighbor interactions and linear dimension 
L = 2000, where a lattice site is chosen at random and the trial spin is set 
to the color of the majority of neighbors. For the pure ferromagnetic one- 
dimensional Ising model, persistence was studied numerically pTj , and for the 
one-dimensional Potts model the exponent was analytically computed [2]. 
For higher dimensions evidence for power law decay was obtained numeri- 
cally [B], while the exponent was analytically obtained [H E] showing that 
P{t) decays algebraically as 

P{t) oct-'^''\ (1) 

where 6{d) is a a non-trivial persistence exponent. For the same ferromag- 
netic Ising model at higher dimensions [6] (c? > 4) and for the two-dimensional 
ferromagnetic Q-states Potts model [S] (Q > 4) it has been shown through 
computer simulations that P{t) is not equal to zero for t ^ oo. This char- 
acteristic is sometimes called "blocking". In this way, if P(oo) > we can 
reformulate the problem by restricting our observation to those spins that 
eventually flip. Hence, we can consider the behavior of the residual persis- 
tence given by 

p{t)=P{t)-Pioo). (2) 

Recently attention has been turned to the persistence problem in systems 
containing disorder [TOl|TT]. Numerical simulations of the zero-temperature 
dynamics of bond diluted (weak dilution [TUI or strong dilution [12]) two- 
dimensional Ising model also reported "blocking" evidence. Howard has 
found evidence of an exponential decay of the persistence with "blocking" 
for the homogeneous tree of degree three (T) with random spin configuration 
at the initial time [16]. In order to investigate other aspects of the disorder 
effect to the persistence problem, Lima et al. [T7] presented new data for 
the zero-temperature dynamics of the two-dimensional ferromagnetic Ising 
model on a random Poissonian lattice. They showed that this system ex- 
hibits "blocking" which means that the persistence does not go to zero when 
t oo. They showed also that p{t) decays exponentially in the long-time 
regime. Recently S. Jain et al. ^3] studied the persistence problem in a 
random bond Ising model of socio-econo dynamics. In this study they found 
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no evidence of the "blocking" effect. Here in this paper we study the persis- 
tence problem in the Potts model with Q states {Q = 3, 4, 2^°) on directed 
and undirected Ba.Taha:Si-AVoeTt(B A) [H] networks and Erdos-Renyi(ER) [15] 
random graphs using the zero-temperature Glauber dynamics. The results 
show that for directed and undirected ER random graphs the persistence de- 
cays exponentially to zero without evidence of "blocking", but for directed 
and undirected BA networks the persistence decays exponentially with and 
without presence of the "blocking". In both cases no algebraic decay of the 
persistence with the time was observed. 
Model and simulation 

We consider the Q-states Potts model on directed and undirected Baiabasi- 
Albert(BA) networks and Erdos-Renyi(ER) random graphs defined by the 
Hamiltonian 

H=-JY: Ss^,s, (3) 

<ij> 

here J > is a ferromagnetic coupling, the sum is over all neighboring 
interactions included in the system, the spin values are in the range Si = 
1, ...,Q, and Sx^y is Kronecker's delta function. We perform simulation over 
networks with = 10^ sites. A random initial configuration of spins is 
obtained and P{t) is calculated over 10^ and 10^ Monte Carlo steps and 
a quenched average is done over 10 different networks for each Monte Carlo 
step. The zero-temperature Glauber dynamics was used in order to check the 
number of spins that never changed their state up to time t. In this dynamics, 
we start with a random initial spin configuration and allow it to be flipped 
at random or following a given logical sequence. The selected spin. Si, is 
flipped or not according to AEi, where AEi is the energy variation of the 
site i. If AEi < the spin Si is flipped with probability one. If AEi = the 
spin Si is flipped with probability 1/2 chosen at random. Finally, if AEi > 
the spin is not flipped. One Monte Carlo step corresponds to application of 
the above rule for all spins of the network. The system configuration is left 
to evolve until a given Monte Carlo step 10^ and 10^. The number n{t) of 
spins that did not change up to this time t is computed for each Monte Carlo 
step for determination of the persistence probability given by [H [TOl [12] 



PM = (4) 

where n{t) is the sites number that did not flip their state up to time t and 
A^ is the number of sites. 
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Results and Conclusion 

In Fig. 1 wc plot P{t) versus t for Potts models with Q = 3,4, 5, 7, 9, 12, 
24,64,128, 256,512,1024,4096,16384 and 2^° states in iV = 10^ sites and 
m = 2 (coordination probability p = m/N for ER random graphs) on undi- 
rected ER random graphs. The persistence probability P{t) seems to decay 
exponentially to zero for all Q values. For directed ER random graphs the 
results are the same. In Figs. 2 and 4 we plot again logP(t) versus t (for 
t from 1 to 500) for Potts model for same Q values {Q = 3, 4, 2'^°) and 
m = 2 (where here m is also the size of the initial core for generating the 
BA networks) values, for directed and undirected BA networks, respectively. 
We observe that for Q = 3, 4, 5, 7, 9,..., 2^° the persistence probability decays 
exponentially to a constant value P{oo) > for directed BA networks. This 
exhibits the evidence of the "blocking" which means that the persistence 
does not go to zero when t oo. Therefore, for undirected BA networks the 
persistence decays exponentially also to a constant value P(oo) for all finite 
Q {Q < 2^^) of the Potts model investigated when t — > oo, again we have 
evidence the "blocking" for these finite Q values, but for Q = 2^° ("Q = oo") 
the persistence probability decays exponentially to zero. In the Fig. 3 and 

5 we plot the same as Fig. 2 and 4, but now for t from 500 to 3 x 10^ for 
Q values {Q = 3, 4, 5, 7, 9,..., 2^° and Q = 3, 4, 5, 7, 9,..., 2^2) for directed and 
undirected BA networks, respectively, showing that (in this observed inter- 
val) there are no evidences for phase transistions from finite values of the 
persistence probability P{t oo) to zero for any value of Q. For directed 
BA networks, P{t) decays to values that correspond to 1, 2 or 3 survivors 
of 10^ sites. For undirected Q — 2^° states the persistence probabihty de- 
cays quickly to zero, showing that there are survivors for Q < 2^2. In Fig. 

6 we plot double-logarithmically the behavior of p(t) versus time for some 
Q = 3, 4, 9, 24 for undirected BA networks. This result showed the exponen- 
tial relaxation of p{t) — > 0. In Fig. 7 we plot the same variables, but now for 
Q = 3,256,512, 1024,4096. In this picture we observe for Q = 3 an expo- 
nential decay to zero, but for Q — 256, 512, 1024, 4096, we see some kind of 
algebraic decline, followed by an exponential decline for p{t) for longer times. 
In Fig. 8 we plot \n{time) versus ln(Q) for Q = 3, 4, 2^^ values and observe 
the time the system needs to reach the "blocking" state on undirected BA 
networks. Initially we observe a growth of the time for values from Q — 3 
until 2^^ states, after that the time decreases for large Q-states values. The 
system needs only a short time to reach the "blocking" state for large Q- 
states values. The linear best fit gives the slope=0.89(3) for growth of the 
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time for values from Q — 3 until 2 states. 

In summary, we have presented new data for the zero-temperature dynam- 
ics of ferromagnetic Potts models on scale free networks and random graphs. 
For the directed and undirected Erdos-Renyi random graphs the system does 
not exhibit evidence of "blocking" independently of the assumed values for 
Q, that is, the persistence probability P{t) go to zero when t oo. In con- 
trast, for directed and undirectedBA networks the system exhibit evidence of 
"blocking" for Q = 3, 4, 5, 7, 9, 2^° values, here observed, when for t — > oo 
the persistence probabihty does not go to zero. For Q large the "blocking" 
there is with only 1, 2, 3 or 4 survivors of the initial total of the 10^ survivors. 
In the undirected case we observe that the persistence probability P{t) de- 
cays exponentially to values different from zero for Q = 3, 4, 5, 7, 9, 2^^ 
when i — >■ oo, that shows the evidence the "blocking": but for Q = 2^° this 
persistence decays to zero, this could be just a finite-size effect since in the 
other "blocking" cases one may have only few unchanged spins. Here wc also 
study the residual persistence probability p{t) for undirected BA networks, 
and believe that for small values Q = 3,4, 9, 24 this decays exponentially to 
zero, but for Q — 256,512, 1024,4096, we see first some kind of algebraic 
decay, followed by exponential dechne to zero. We also observe the increase 
of the time to reach "blocking" with increasing Q, that initially is linear in 
the log- log scale of the Q = 3, 4, 7, 2^"^ and afterwards the time decays 
dramatically to short times for Q = 2^^, 2^°, 2^^. 
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Potts model on undirected Erdos-Reny graphs 
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Figure 1: P{t) (on logarithmic scale) versus time for various Q and m = 2 
on undirected Erdos-Renyi(ER) random graphs. Q = 3,4, 5, 2^° from top 
to bottom. 
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Potts model on directed Barabasi- Albert networks 

N=10^ sites, Q=3, 9, 12, 24, 64, 128, 512, 2^°, 2^^ , 2^*, 2^^, 2^^ 2^° & 2^° from top to bottom 
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Figure 2: P{t) (on logarithmic scale) versus time for various Q and m = 2, 
on directed Barabasi- Albert (BA) networks. Q = 3, 4, 5, 2^° from top to 
bottom. 
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Potts model on directed Barabasi- Albert networks 

N=10^ sites, Q=3, 9, 12, 24, 64, 128, 512, 2^°, 2^^ , 2^*, 2^^, 2^^ 2^° & 2^° from top to bottom 
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Figure 3: P{t) (on logarithmic scale) versus time for various Q and m = 2, 
on directed Barabasi- Albert (BA) networks. Q = 3, 4, 5, 2^° from top to 
bottom. The time varies from 500 to 300000. 



9 



Potts model on undirected Barabasi- Albert networks 

N=10^ sites, Q=3, 9, 12, 24, 64, 128, 512, 2^°, 2^^ , IL^, IL^ , 2^^ 2^° & 2^° from top to bottom 
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Figure 4: Pit) (on logarithmic scale) versus time for various Q and m = 2 
on undirected Barabasi- Albert (BA) networks. Q = 3, 4, 5, 2^° from top to 
bottom. 
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Potts model on undirected Barabasi- Albert networks 

N=10'^ sites, Q=3, 9, 12, 24, 64, 128, 512, 2^°, 2^^ , 2^^ 2*^ 2^^ & 2^" from top to bottom 
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Figure 5: P{t) (on logarithmic scale) versus time for various Q and m — 2 
on uncJzrecierf Barabasi- Albert (BA) networks. Q — 3, 4, 5, 2^^ (because 
Q = 2^° is zero in this interval) from top to bottom. The time of 500 to 
300000. 
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Potts model on undirected Barabasi- Albert networks 

N=10^ sites, Q=3, 4, 9 & 24 from bottom to top 
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Figure 6: p{t) versus time (double-logarithmic scales) for various Q = 
3,4,9,24 (from bottom to top) and m = 2 on undirected Barabasi- 
Albert(BA) networks. 
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Potts model on undirected Barabasi- Albert networks 

N=10^ sites, Q=3, 256, 512, 1024 & 4096 from top to bottom 
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Figure 7: p{t) versus time (double-logarithmic scales) for Q = 3 and other 
various values of Q (256,512, 1024,4096) (from top to bottom) and m = 2 
on undirected Barabasi-Albert(BA) networks. 
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Potts model on undirected Barabasi- Albert networks 

16 



XT in^ a /I n O/i o** o'^ o'" o'^ o'^ 16 18 19 20 22 

N=10 Sites, Q= 3,4, 9, 24, 2, 2, 2, 2 ,2 ,2 ,2 ,2 ,2 ,2 &2 



14 

12 
10 
5, 8 
6 
4 
2 




' 1 ' 


1 1 1 1 1 1 1 1 1 1 1 1 1 




o ° ° 




o 




o 




o 




o 




o 


o 




o 




cP 






o o o 


1 


1 . 1 . 1 . 1 . 1 . 1 . 



2 4 6 8 10 12 14 16 18 

Ln(Q) 



Figure 8: In(time) versus ln{Q) for various Q = 3, 4, 5, 2^^ and m = 2 on 
undirected Barabasi-Albert(BA) networks, where fAme is the number of time 
steps before "blocking" occurs. Here we observe that time increase when Q 
increases for Q — 3, 4, 5, 2^^. The best linear fit gives a slope=0.89(3). 
Afterwards the time decays dramatically to short times for Q — 2^^, 2"^^, 2^^. 
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